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Abstract 

Correlations between non-identical particles at small relative velocity probe 
asymmetries in the average space-time emission points at freeze-out. The 
origin of such asymmetries may be from long-lived resonances, bulk collec- 
tive effects, or differences in the freeze-out scenario for the different particle 
species. STAR has extracted pion-proton correlation functions from a dataset 
of Au-fAu collisions at ^Js^n = 200 GeV . We present correlation functions 
in the spherical harmonic decomposition representation, for different central- 
ities and for different combinations of pions and (anti-)protons. 



1 Introduction 

Femtoscopy is the primary technique that provides possibility of measuring 
space-time quantities of the order of 10~^^m and lO^^^s. It allows for inves- 
tigation of evolution of nuclear matter created in heavy ion collisions. Two 
charged particles with small relative velocity emitted from source interact 
through strong and Coulomb interactions called final state interactions (FSI) 
which give correlation between them. Correlations of non-identical particles 
are sensitive to space-time asymmetry in emission of particles |lj |2j- By mea- 
suring them we can investigate which sort of particles are emitted earlier or 
later and what is the relative shift between average emission points. 

Information about size and space-time asymmetry of the source can be 
extracted from correlation functions calculated with respect to the sign of 
the out (along the velocity of the pair in transverse plane), long (along the 
beam axis) and side (in transverse plane, perpendicular to the other two) 
components of the vector k* (momentum of first particle in Pair Rest Frame). 
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Such preliminary pion-proton correlation functions measured by STAR were 
presented in [3]. However such information is most efficiently encoded in a 
spherical harmonics representation of the correlation function [1]. In case of 
limited momentum acceptance of the detector or big difference in mass be- 
tween particles, when some orientation in space of the k* vector are not regis- 
tered by the detector, calculation of correlation function directly in spherical 
harmonics is the most effective method [5]. Functions calculated directly in 
spherical harmonics are less sensitive to momentum acceptance limitations 
than functions calculated from 3-dimensional correlation function. 

2 Experimental procedure 

This work presents analysis of data taken in 2004 by the STAR experi- 
ment. Au+Au collisions at ^/SNN = 200 GeV are analyzed in 3 different 
centrality bins: 0-10% - central events, 10-30% - intermediate events and 
30-50% - mid-central events. The particle identification is done in Time 
Projection Chamber (TPC). Both pions and protons are selected based on 
Ncr = [—3.0, 3.0] value which describes distance from mean of the Gaussian 
parametrization of the dE/dx curve. Particles are selected from the mid- 
rapidity region \y\ < 0.7. Transverse momentum of the pion is set within the 
range of pt = [0-1, 0.6] GeV jc and for proton = [0.4, 1.25] GeV jc. 

Pion-proton pairs are constructed from close velocity particles. Due to 
significant mass difference between pions and protons, pion with momentum 
0.18 GeVjc requires proton with momentum 1.2 GeVjc. This shows that 
this analysis reaches the limits of the particle identification capabilities of the 
STAR TPC detector. At the level of pair selection, pairs of tracks which share 
more than 10% of the hits in the detector are removed from the analysis (the 
so called hit merging effect). We also eliminate correlated electron-positron 
pairs which come from gamma conversion (based on topological cut) and non 
pion-proton pairs (based on probability). 

Background is constructed only from events with similar characteristics 
thus events are divided into bins with respect to z-vertex position - 15 bins, 
event multiplicity - 6 bins and mean transverse momentum of the registered 
particles - 3 bins. 

To obtain experimental correlation function we construct distribution of 
the correlated pairs (coming from the same event) and distribution of the 
uncorrelated pairs (coming from different events). Distributions are directly 
binned in spherical harmonics. Detailed explanation of the formalism can be 
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found in [5]. 

In our analysis CQ{k*) and ReC\{k*) components are analyzed. Cq rep- 
resents overall size of the system. ReC\ component represents space-time 
asymmetry observed in the out direction. 

Constructed functions are corrected on purity which is defined as a prod- 
uct of PID probability of the particles in the pair and fraction of primary 
pairs. 



3 Results 

Figures [1] and [2] show Cq and C\ components of the correlation function for 
the — p and 7i~ — p pairs correspondingly calculated for Au+Au colli- 
sions at i/J/wv = 200 GeV. Red dots show central data, blue squares are 
for intermediate events and green triangles are for mid-central events. Peak 
around A;* = 0.1 GeV/c reflects correlated products of the (anti)lambda de- 
cay. Lambda peaks do not overlap with femtoscopic correlation effect. Cq 
functions show positive correlation as the interaction between unlike sign par- 
ticles is repulsive. ReCl functions show small deviation from zero for small 
values of k* which suggests that pions and protons are not emitted from the 
same regions of the source or are not emitted at the same time. Figure [3] 
presents functions for like sign pairs where correlation effect is negative as 
a Coulomb interaction is attractive. C} components for like sign pairs also 
show asymmetry. 
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Figure 1: — p correlation functions. Left panel - Cq component, right 
panel - ReC} component. Red points - central 0-10% events, blue points - 
intermediate 10-30% events and blue points - 30-50% mid-central events. 



In the fit procedure 3-dimensional Gaussian profile of the source is as- 
sumed. Set of correlation functions with various parameters of the source 
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Figure 2: 71^ — p correlation functions. Left panel - Cq component, right 
panel - ReCl component. Red points - central 0-10% events, blue points - 
intermediate 10-30% events and blue points - 30-50% mid-central events. 



profile is calculated Then around the best value an ellipse of co vari- 
ance is calculated. Figure H] shows results for experimental data as well as 
results obtained for Therminator-|-Lhyquid[7j and UrQMD v3.3[8] data. Ex- 
perimental results obtained for unlike-sign pairs show centrality dependence 
of the size {(Tout) and asymmetry (fiout) of the source. Model study is done for 
both like-sign and unlike-sign pairs. In UrQMD data observed asymmetry 
seems to be constant for investigated centralities. 

Experimental procedure allow only for measuring total asymmetry, more 
detailed study with models gives more detailed information. Figure [5] shows 
distribution of the emission time difference versus separation of the emis- 
sion points calculated in Therminator and UrQMD. Positive values of ryj^e 
corresponds to pion emitted later than proton. Such scenario dominates in 
Therminator data. Negative values of Tout means that emission point of pro- 
ton is shifted outward (to the edge of the source) with respect to the emission 
point of pion. 



4 Summary 

Analysis presented here shows that in system created in Au-|-Au collisions 
at ^/sNN = 200 GeV average emission points of pions and protons are not 
the same. Observed asymmetry depends on centrality and is correlated with 
the size of the observed pion-proton system. Technique of data analysis al- 
lows only for measuring total asymmetry but study done with the models 
gives possibility to estimate what is the contribution of emission time dif- 
ference and pair separation to the total asymmetry. Such study done with 
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Figure 3: 7r+ —p (up) and 7r~ — p (bottom) correlation functions. Red points - 
central events, blue points - intermediate events and blue points - mid-central 
events. 



Therminator suggests that protons are mainly emitted earlier than pions. In 
UrQMD such phenomena is not observed. Both models show that in average 
protons are shifted to the edge of the source. 
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Figure 4: Ellipses of covariances show the 3D Gaussian fit results of the ex- 
perimental (left) and theoretical (right) correlation functions. Experimental 
results are for two orientations of the magnetic field in the detector (dash-dot 
line). Continuous lines show averaged results. 
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Figure 5: Distribution of the emission time difference versus separation of 
the emission points in Therminator (left) and UrQMD (right). 
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